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In this work we study the effects of collapse and revival as well as Zitterbewegung (ZB) phe- 
nomenon, for the relativistic electron wave packets, which are a superposition of the states with 
quantum numbers sharply peaked around some Landau level no of the order of few tens. The prob- 
ability densities as well as average velocities of the packet center and the average spin components 
were calculated analytically and visualized. Our computations demonstrate that due to dephasing 
of the states for times larger than the cyclotron period the initial wave packet (which includes the 
states with the positive energy only) loses the spatial localization so that the evolution can no longer 
be described classically. However, at the half-revival time t — Tr/2 its reshaping takes place firstly. 
The behavior of the wave packet containing the states of both energy bands (with E n > and 
E n < 0) is more complicated. At short times of a few classical periods such packet splits into two 
parts which rotate with cyclotron frequency in the opposite directions and meet each other every 
one-half of the cyclotron period. At these moments their wave functions have significant overlap 
that leads to ZB. At the time of fractional revival each of two sub-packets is decomposed into few 
packet-fractions. However, at t — Tr each of the two sub-packets (with positive or negative en- 
ergy) restores at various points of the cyclotron orbit, that makes it impossible reshaping of initial 
wave packet entirely unlike the wave packet which consists of states with energies E n > only. 
*^ ' Obtained results can be useful for the description of electromagnetic radiation and absorption in 

relativistic plasma on astrophysics objects, where super high magnetic field has the value of the 
order 10 8 — 10 9 T, as well as for interpretation of experiments with trapped ions. 
Ch 
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INTRODUCTION 



The space-time evolution of quantum system with discrete but not equidistant energy spectrum can be quite 
complex. It was showni that for the initial state consisting of the eigenstates with quantum numbers n sharply 
peaked around some level no ^> 1 the coefficients of Taylor expansion of spectrum E n around that value define the 
£C) ■ different time scales of wave packet periodic evolution. 

At the beginning a localized wave packet moves periodically along a classical trajectory with the period T c \ = 
27rh/E' nQ . However, over time, the nonlinear terms in (n — no) in the expansion of the spectrum E n become significant, 
resulting in a spreading of the packet. The next stage of a long-term evolution of semiclassical wave packets for a large 
class of quantum systems is universal as was shown in Ref. 1 . Thus, at the revival time Tr = A-kTi/ E'^ o the additional 
phases of the individual components of the wave packet due to quadratic terms in the expansion are exactly multiplied 
by 27r which means a full reshaping of the wave packet. At the times of fractional revivals t = mTn/n, where m and 
^ n are mutually prime integers, the phase differences between subsets of the component eigenfunctions are stationary, 

and the wave packet breaks up into sets of sub-packets. After a few revivals one obtains a complicated distribution 
of probability density around the classical trajectory. 

During some past decades the phenomena of collapse and revival have been investigated theoretically and experimen- 
tally in different quantum systems. In particular, revivals and fractional revivals have been investigated theoretically 
for Ridberg atoms, molecules and in low dimensional quantum structures. Firstly, the detailed calculation of a short 
laser pulse exitation of the high-energies Rydberg states and the total power radiated by the atom was performed in 
Ref. 2 . The space-time dynamics of the hydrogen atom was considered in Ref. 3 by constructing a minimum- uncertainty 
wave packet that travels along a kepler orbit. It was shown that both classical and quantum properties associated 
with the discreteness and non-equidistance of the atomic spectrum occur inevitably in the long-term evolution of 
wave packets. Romera and Santos 4 studied the cyclotron dynamics of electrons in a monolayer graphene, where low 
energy excitations are massless Dirac fermions. They showed that when electrons are described by wave packets with 
a Gaussian population of only positive energy Landau levels, the presence of the magnetic field induces revivals of the 
electron currents, besides the classical cyclotron motion. When the population comprises both positive and negative 
energy Landau levels, revivals of the electric charge current manifest simultaneously with Zitterbewegung (ZB) and 
the classical cyclotron motion. In Ref. 5 the quasiclassical evolution as well as the phenomenon of collapse and revival 
of a wave packets in graphene quantum dots in a perpendicular magnetic field was studied. 



The Dirac equation also predicts an unexpected features of relativistic motion. Bermudes et al. 6 considered the 
evolution of relativistic wave packets in a magnetic field and predicted three regimes: macroscopic, microscopic and 
mesoscopic (intermediate) when the electron energy level comprises several tens of Landau levels. The authors of 
Ref£ applied an exact mapping between relativistic model and a combination of Jaynes-Cummings and anti-Jaynes- 
Cummings interactions 7 widely used by the quantum optics community. In the work 6 the relativistic version of 
"Schrodinger cat" states, or "Dirac cat" states, was built for relativistic Landau levels when an external magnetic 
field couples the states of a relativistic spin 1/2 charged particle. 

Rusin and Zawadzki^ considered the effect of trembling motion of wave packet center or ZB for relativistic electrons 
moving in a vacuum in the presence of an external magnetic field. As in the original work by Schrodinger 9 they used 
one-electron approximation and found that in this case the effect of ZB is very small. The authors of Ref. 8 noted, 
that in accordance with results obtained by Krekora 10 , the fully occupied electron negative energy states, or Dirac 
vacuum, may prevent the interference and ZB to occur. For this reason the Zitterbewegung is merely a mathematical 
property of the one-electron Dirac equation and can not be observed for real electrons. So, the calculations in Ref. 8 
were made for 2+1 and 3+1 Dirac equations for parameters, which correspond to the motion of trapped ions in recent 
experiments by Gerritsma and coworkers^. 

Gea-Banacloche evaluated the asymptotic description for the evolution of two-level atom interacting with a quan- 
tized field in initially coherent stated. The collapse and revival was investigated in this system. It was stated that 
collapse appears to be associated with a measurement of the initial state of the atom with the field. At the half-revival 
time the pure state of the field, which is a mesoscopic superposition of coherent states, is referred to as a "Schrodinger 
cat" . 

The phenomena of collapse and revival was observed experimentally in different nonlinear quantum models. Among 
them, the two-level atom trapped in a high-Q microwave cavity 13 , Rydberg wave packets in atorn^ and molecular vi- 
brational states. 15 Interestingly, methods for laser isotope separation which take advantage of the revival phenomenon, 
were proposed and realized in Ref. 16 . In experiment Monroe et all 17 a superposition of two different coherent states 
was created for an ion oscillating in a harmonic potential. 

The quantum dynamics of free relativistic particles represented by three-dimensional Gaussian wave packets with 
different initial spin polarizations was investigated and visualized in Ref. 18 . The influence of the symmetry of the 
initial wave packet on its kinematic characteristics such as average electron velocity, the direction of trembling motion 
as well as the distribution of spin densities was analyzed also. The effects of splitting and ZB of electron wave packets 
in the semiconductor quantum well under the influence of the Rashba spin orbit coupling and also in monolayer 
graphene were considered in Refs.— »2£. 

In this work we study the long-term cyclotron dynamics of relativistic electrons, which obey the Dirac equation. For 
this purpose we use the inertial frame in which the component of momentum parallel to the direction of the magnetic 
field equals to zero. Such choice allows to represent the most of results in the clear analytical form. To describe the 
behavior of the relativistic particle we consider the superposition of the eigenstates of the Dirac Hamiltonian which 
is centered around some quantum number no ^> 1 . The bispinor structure of electron wave function and existence of 
two energy bands lead to a quite complex dynamics of wave packets. The space- time evolution of electron probability 
densities, as well as spin densities was investigated analytically and numerically and after that visualized. 

The paper is organized as follows. In Sec. II the relativistic electron eigenstates and eigenvalues of Dirac Hamiltonian 
in a magnetic field are introduced and the analogy with Jaynes-Cummings model is discussed. In Sec. Ill we consider 
the dynamics of electron wave packets when the population of Landau levels comprises positive energy states only. 
The structures, which form at fractional revival times are visualized and physical effects of appearance of these 
structures are discussed. The detailed description of spin precession is presented. In Sec. IV we consider the evolution 
of mesoscopic wave packet, which is a superposition of the states with positive and negative energies. We discuss some 
peculiarities accompanying the revival phenomenon for such wave packet, the effect of ZB, in particular. The results 
obtained in this Section can be useful for observation of collapse and revival in experiments with trapped ionsi^. 
Finally, in Sec. V, we conclude with the discussion of results. Some auxiliary mathematical results are presented in 
Appendices A, B, and C. 

II. MODEL AND APPROACH 

We consider the cyclotron dynamics of relativistic electron in a magnetic field B = (0,0,5). Using the Landau 
gauge, where A = {—B ■ y, 0, 0), the Hamiltonian can be expressed as follows: 

H = ca x (p x - eBy/c) + ca y p y + 

+ca z p z + mc 2 /3, (1) 



where c is the light velocity, m is the electron mass, (— e) is the electron charge and c^ and f3 are Dirac matrices 

«i=(° o)' P= (o -/)' i = 1 ' 2 ' 3 ' ( 2 ) 

cr^- are Pauli matrices. The Dirac equation admits a solution of the type 



( ip x x+ip z z \ 

\ n ) 



exp 
*(r) = \ nh '- U{y). (3) 

For simplicity, we assume in Eqs.(l) and (3), p z = (we always may choose the inertial frame, in which this condition 
takes place). Then the equation that we consider is 

HU = EU, (4) 

where the effective Hamiltonian H 

ft = ( mc2 cPyVy ~ eB (v ~ Vc)crx \ ,^ 

\cpydy -eB(y-y c )a x -mc 2 y' i; 

and y c = cp x /eB. It is easy to check that the operator 

commutes with Hamiltonian (5) and consequently we can construct a function which is an eigenfunction of H and K 
simultaneously. As follows from Eq.(6), the eigenstates of K with eigenvalues A& = ±1 can be written as 

</> Afc =i = (/| t>, x\ 4>) T = (/, 0, 0, X f, (7a) 

V> Afc= -i = (/| i>, X \ t>) T = (0, /, X, Of, (76) 

where / = /(r, £), \ — x( r ?^) are arbitrary functions of r and t. Spin-up | \> and spin-down | l> spinors are given 
by 

it>=(;), u>=(;)- w 

We see that for a given value A& two components of the Dirac spinors (7a), (7b) are vanished: ^2 = ^3 = for 
Afc = 1 and ipi = ^4 = for A& = — 1. Thus, depending on the sign of A&, the Hilbert space of solutions of the 
Dirac equation (4) with Hamiltonian Eq.(5) can be divided into two invariant subspaces. Correspondingly, the Dirac 
Hamiltonian is decomposed into two terms. One of them couples the components ipi and ^4, and this part is identical 
to a Jaynes-Cummings model, describing the interaction of two-level atom with a quantized single-mode field^ The 
remaining term which couples the components ip2 and ips is identical to the anti- Jaynes-Cummings interaction^ This 
suggests that we can deal with two-component wave function ip(y) = (f(y),x(y)) T - For A& = 1 one obtains from 
Eqs.(4), (5) and (7a) 

mc 2 -V2he~c~Ba\ ( f(y)\ _ „ ( f(y)\ (q , 

-V2he~c~Ba^ -mc 2 J \ X (v) )~ \ x(v) J ' l j 

Here we introduce the raising and lowering operators for the harmonic oscillator 



» -Uf?* -*-£=)• 



'-Mf?o>-^&;)- 



where uj = eB/mc is the cyclotron frequency for a non-relativistic electron. The eigenvalues of Eq.(9) are: 



E s ^ n = s • e n = s\fm 2 c A + 2TteBcn (11) 

with n = 1, 2, 3, . . . for s = 1 and ra = 0, 1, 2, 3, . . . for 5 = — 1. Correspondent eigenstates are given by 

^.-...w-^fi':^)' (12) 

where (j) n (y — y c ) is linear oscillator wave function, and the coefficients d n and b n are determined by 



Two-component wave function for A& = — 1 obeys the Dirac equation which is similar to Eq.(9): 



mc 



-V2h^Ba^\ ff(y)\ = „ ( f(y) 



e x • (is) 



-\/2hecBa -mc 2 J \x(y) J \x(y) J ' 

As before, the energies E s ^ n are determined by Eq.(ll), but now n = 0, 1, 2, 3, . . . for s = 1 and n = 1, 2, 3, . . . for 
s = — 1. The associated eigenstates are 

*■-■--<»>-(-&&))• (i6 > 

Note that in the considered case p z = the four-component wave function [/(?/) in Eq.(3) is connected with two- 
component wave function ^\ jS ,n(y) (see Eqs.(12), (13), (16) and (17)) by relations 

U Xk = hs ,n(y)=( l l > 1^)^=1,^(2/), (18) 

U Xk = - 1 , s ,n(y)=( li Q > ^^--l.sAv), (19) 

Thus, in the considered case (p^ = 0) the general solution can be written as 

i/;(r,t)= dp<p p (x) J2 C n,x k ^ U ^,s,n(y) e ^(- is£ nt/h), (20) 

n,s,Afc 

where cp p (x) — l/V27rhex.-p(ipx/h), and coefficients C* x (p) are to be determined by the initial wave function. 

III. EVOLUTION OF WAVE PACKET FORMED BY THE POSITIVE ENERGY EIGENSTATES 

A. Time dependence of the electron probability density 

To study the complex dynamics of a real relativistic electron placed in the magnetic field, we construct the initial 
mesoscopic wave packet as a linear combination of the positive energy states only. Besides, we consider the special 



form of the initial wave function which represents the coherent state of the non-relativistic electron in the magnetic 
field. So, let the coefficients C^ x (p) in Eq.(20) are determined by the expressions: 



s = l ,„\ _ Pg(P) C n+l n8=1 ,s _ OLg(p)Cn 



Cn>0,\ k = -l\P) — / 9 ; n ^ i Cn>l,A fe = l(P) 

ys= — l 



^TW ™^=— ^TW 

cta; 1 (p) = ' ( 21 ) 



where 

g(p) 



\ ft~7^ exp v 2^ pa ^ n ~ qa ^)' ^ 



= exp(-(gq) 2 /4) • (-gq)"- 1 



a = ^/hc/eB is the magnetic length, the parameter ga ~ y^no characterizes the radius of relativistic orbit. For 
simplicity we assume a = a*, f3 = f3* . Then from Eq.(20) we have 



^fo^ = / 9 o9 / d P ( P P ( X )9(P) X 

^OL 1 + (3 2 J 

y^ ( d n (acn0 n _i(2/)| t> +/3c n+ i0 n (2/)| |>) \ -ie n t/h , 
Z-j I -6 n (acn0 n (2/)| |> +/3c n+ i0 n _i(?/)| t>) / 



(24) 



Note that in the considered case all components of the wave function differ from zero. Thus, for describing the wave 
packet dynamics we should use the full Hamiltonian, Eq.(l). To analyze the motion of the packet center we first of 
all have to find the average value of the velocity operator V{ = cai. Using Eq.(2) and Eqs.(22), (23), (24) we arrive 
after integration over x and y at the following equation: 



, N ce~^ I 2 v^ M 2n+1 / (p n+1 - 1 



a 2 hfn+2 + l c i(y n + 2-V n + l)T + 02 r£j!^tl e K<Pn + l-<Pn)T j _ ^5) 



^n+2 V ^n 

Here and below y? n is the electron energy (in units mc 2 ) 



Vn = Vl + 2n(A/a) 2 , (26) 

and r = ct/A, where A = fi/mc is the Compton wave length. It can be shown that the dominant contribution to 
the sum in Eq.(25) comes from the interval in the neighborhood of n « no = (qa) 2 /2. As follows from Eq.(26) for 
n(A/a) 2 <C 1, ip n ~ 1 + n(A/a) 2 and Eq.(25) describes the cyclotron motion of a non-relativistic electron 

fog 

^cU) = COSCjt, 

m 

vZ(t) = — smut. (27) 

y m 

If no is large enough, the quantum interference of the wave packet components leads to the different types of 
periodicity. In Fig.l we plot the dependence oiv^(t) for a relativistic wave packet with qa = 5, a = (3, placed in the 
magnetic field B ~ 4.5 • 10 T T. (Note that this magnetic field is 1.5 times greater than the maximum value achieved 
currently in laboratory experiments). In this case (f n +i — Vn ~ ^p' no + (^ /2)(2n — 2no + 1). So, in the beginning 
the average velocity of the wave packet center oscillates with classical cyclotron motion period 

2irh 2iTip no eBc 

J-Cl = o— f— = 5 ^C = • [*0) 
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FIG. 1: The dependence of the average velocity v x (t) versus time for the initial wave packet, Eq.(24), with a = /3 and the 
parameters A/a = 0.1, qa = 5 (that correspond to the value n = no w 13 and dispersion An ~ ^/no « 3.5). 



At large times t > To these oscillations fade away but then appear again. The damping time Tp can be estimated 
from Eq.(A4). For |^" |r <C 1 the amplitude of the oscillations is proportional to exp( — (qa(p , / lQ r/2) 2 ) = exp(— t 2 /T^), 
where 



T D 



2A 



2% 



<iac\<Pn \ 



qa>K 



(29) 



At times Tjj < t < Tr/2 the average velocity components as well as average coordinates of the packet center are 
negligible. At t ~ Tr/2 (see Fig.l) the velocity oscillations will occur again. Here Tr is the revival time at which the 
wave packet is fully reformed. 1 At tk = kTj$ where k is integer, the nonlinear (quadratic) term in the expansion of ip n 
has no effect, so that 



Tr 



4ttA 



A-kTi 



(30) 



Note that according to Averbukh and Perel'man 1 the moment Tr/2 is the first reconstruction time of the wave packet 
(see more detailed discussion below). For the considered wave packet Eqs.(29) and (30) give Tr/Tjj ~ 30 that is in a 
good agreement with numerical calculations (Fig.l). 

Now we describe some peculiarities of the space-time dynamics of the electron wave packet (24). Performing the 
integration over p in Eq.(24), one has 



^(r,t) 



E 



7 



n-l 



7 



(1 - S nt o)e %e (adn, 



pb n qa ab n p _ ie \ T 



2n 



2n 



^(0,1,0,0^ 



-idn—icp n ct/X 



The details of the calculation are presented in the Appendix B. Here we use the notations 

.2 



X 

a 



. y — qa* 
psmt/, = pcostf, 7 



qap 
"~2"' 



(31) 



(32) 



M(p,0) = 



exp l 



p 2 -\-{qa) 2 , ipsin 9(pcos 9-\-2qa) \ 
4 ' 2 J 



a\/2ir 



(33) 
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FIG. 2: (Color online) The probability densities for the wave packet, Eq.(24), with a = f3 and the parameters A/a = 0.1, 
qa = 5 at times: (a) t = 0, (b) t w T D , (c) t « Tr/4, (d) £ w T^/2, (e) t»T fl . 



In Fig. 2 we plot the time evolution of the probability density \?p(r, t)\ 2 calculated with the use of Eqs.(31), (32) and 
(33). 

At the initial stage t <C Tjj the motion of the wave packet is determined by the linear term of the expansion (A2). 
In this approximation 



i/>(r,t) =i/; c i(r,t), 



(34) 



where i/j c i(r,i) is determined by Eq.(B4). As was noted above in our case the value 2no(A/a) 2 ~ (q^) 2 = 0.25 so 
that ijj c i(Y,t) differs slightly from the coherent wave function of non-relativistic electron (see Eq.(B4)). Therefore, 
the shape of the wave packet \ip c i(r,t)\ 2 remains unchanged during its cyclotron motion, Eq.(B6), with the classical 
period T c \ of the order of attoseconds (10 -18 s.). 

However, after a few orbits at t ~ Tjj (Fig. 2b) the quadratic term in the expansion (A2) becomes significant, that 



leads to the dephasing of different terms in the superposition Eq.(31). As a result, the wave packet is decomposed 
in a host of sub-packets. As was shown by Averbukh and Perel'man 1 at time t ~ tuTr/h, where m and n are 
integer and m/n is an irreducible fraction, the number of such sub-packets is equal to N = n(3 — (— l) n )/4. The first 
reconstruction of wave packet takes place around time t\ = Tr/2X where the revival time Tr is determined by Eq.(30) 
and Tr « 2.3 • 10 -16 s in our case. At that time t\ the n— th term in the superposition (31) acquires the phase factor 
fk = e l7rk = e 27rfc , k = n — no, that leads to the result for the wave packet function (see Fig. 2d) 

^(r,ti)=^(r,ti + -^). (35) 

It should be noted in accordance with Eq.(B7) that the position of the reconstructed packet at the moments Tr/2 
and Tr does not coincide with initial position of the packet at t = 0. At early times, for example, at t^ = Tr/4, the 
original wave packet splits into two subpackets (the fractional revival). In fact, the additional phase factor in this 
case fk = e l7rk I 2 is a periodic function fk = //c+2, so one can use the ordinary Fourier expansion 1 

fk=Po+Pie iirk ,k = 0,l, (36) 

from which 

e i7r/4 e -i7r/4 

P0 = ^ T ,p, = ^ 7r . (37) 

Then, by substituting Eqs.(36) and (37) into Eq.(33) we obtain 

ITT /4 — 27T/4 

^(r,t 2 = Tr/A) = l_^(r,t 2 ) + i-_^(r,fa +Tcz/2). (38) 

Note that according to Eqs.(B7) and (38), the angular distance between these two packets is 180° (Fig. 2c). Similar 
structure consisting of N sub-packets evenly distributed along the classical orbit arises at the other fractional revival 
times. It should be noted, however, that the effects of fractional revivals are not manifested in the dependence of the 
mean velocity on the time (see Fig.l). The shape of the wave packet at revival time Tr (Fig.2e) differs from the initial 
one (Fig. 2a). The spreading of the packet at t = Tr is connected with the next (cubic) term in the Taylor series. 

Features of quantum dynamics of the relativistic wave packets described above should be seen in the nature of 
electromagnetic radiation emitted by the moving electrons. The radiation field as well as its intensity are determined 
by the multipole moments of system. Thus, at the initial stage of the evolution of wave packet (t < To), the dipole 
radiation with a frequency uo = uj c = ^^ is dominated. Its intensity (in the classical approach) is proportional to the 

square of the second time derivative of the average dipole moment d(t) = — ev(t), where the components of average 
velocity are determined by Eq. (25). At times Tjj < t < Tr/2 the mean velocity of the packet is negligible (Fig.l). 
During this time electromagnetic emission is determined by both the time-dependent components of the quadrupole 
moment tensor D a ^ = f |\I/(r, t)\ 2 (3x a xp — Y 2 5 a ^)dY, and the magnetic moment m = ^ f \£(r, £)t[r, v]\£(r, t)dr, 
where the wave function \£(r, t) is determined be the Eqs. (31)-(33). In particular, at the time moment t = Tr/4, 
when the initial wave packet splits into two sub-packets, the quadrupole radiation at the doubled frequency 2u c is 

dominated, that is confirmed by numerical calculations of radiation intensity along n direction: dl /d£t = 1A l^ c s [D, n] 2 . 
In the time interval At ~ 2T]j near the moment t ~ Tr/2 (see Fig.l) the dipole radiation begins to dominate again. 
Then at large times the multipole and magnetic-dipole terms become significant, and so on. 

B. Spin dynamics 

Below we describe and visualize the spin dynamics of the electron wave packet, Eq.(24), in the presence of a magnetic 
field B = (0,0,5). It is well-known that in non-relativistic quantum mechanics the electron spin precesses in the 
(x, y)-plane with cyclotron frequency uo = eB/mc, so that the ^-component remains constant. The spin operators S 
for the Dirac particle 

5, = 5^ = 5fn °.V * = 1,2,3. (39) 



2 l 2\0 a, 

do not commute with the Dirac Hamiltonian, Eq.(l). It follows that besides the usual spin precession which is more 
complicated in the relativistic case, the ^-component changes with time in general. However, the average value of 



S y (t)/(ft/2) 



Tr/2 




FIG. 3: Time dependence S x (t) for electron wave packet, Eq.(24), with the parameters A/a = 0.1, qa = 5. 



component S z is conserved for the wave packet consisting of the states with E n > or E n < only. It directly follows 
from the fact the commutator [S z , H] is equal to zero in the subspace of the eigenvectors of the Dirac Hamiltonian 
for the same sign of energy. The original wave packet (Eq.(24)) satisfies this condition. Thus, we are interested to 
find the time-dependent average values of S x (t) and S y (t) 



Sx,y{t) 



dY^ ] (Y,t)Yi x ^(Y,t), 



(40) 



where ip(r,t) is given by Eq.(24). It is convenient to represent the result of calculations in the complex form 



£„.(*) + »E„(t) 



-/3 ; 



n=0 




n(<p n - i)(<p n - 



+ 1 



i: 



(n + l)ip n ip n 



+1 



ict(ip n+1 -tp n )/\ 



(41) 



The dependence S x (t) represented in Fig. 3 is similar to the behavior of the average velocity component Vx~{t) (see 
Fig.l). Thus, we arrive at the conclusion that the usual precession continues during several periods T c \ only. At 
time t < Td the spin-rotation stops and arises again later at t <J Tr/2. Note that the spin precession of a moving 
electron in a uniform magnetic field is closely connected to the conservation of helicity E^ • (p + eA/c)j. Notice also 
that, as it follows from Eq.(41), the average values of spin components E^ = E y = if the electron wave packet is a 
superposition of the states with A& = 1 (or A& = —1) only (see Sec. I). 

Now we consider the spin densities which we denote as < S Xi y(r,i) >= | < E^(r,£) >, where 



< Z x ,y(r,t) >= ^ f M)E^M). 



(42) 



Using the expression for the wave function, Eq.(24), we find from Eq.(42) < E x (r, t) > and < E y (r,t) > as the real 
and imaginary parts of 



a Q e -ttqa) 2 +p 2 )/2 

<Z x (p,0,t)>+i<Z y (p,0,t)>= Z_ 2 ,^ 2 , p2 , Yl 



27ra 2 (a 2 +(3 2 ) 



m,n=0 



Hfgg/g) 
n\m\ 



n-\-m 
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FIG. 4: (Color online) The distribution of spin density < T, x (p,9,t) > for the wave packet, Eq.(24), with the parameters 
A/a = 0.1, qa = 5 at times: (a) t = 0, (b) t » T R /4, (c) £ » T H /2, (d) £ w T R . 



/ (^m+l +l)(^n + 1) / n((^ n - l)(y? m +l ~ jQ \ 

exp(zct((/p m+ i - ip n )/\ + i<9(ra - n)). (43) 

Fig. 4 and Fig. 5 illustrate the distributions of < Y, x (p,0,t) > and < T iy (p^0^t) > at the same moments as an Fig. 2, 
i.e. t = 0, Tr/4, Tr/2, Tr. Just as for electron probability density at £ < Tr in a main approximation the space-time 
evolution of the spin density is determined by the wave function ip c l( r : t), Eq.(B4), which corresponds to the "classical" 
motion of the electron. So that let us define the functions (a = j3) 



<E?(p,M)>=exp(- 



p 2 + (qa) 2 + 2pqa cos(6» + 2nt/T c ) • 



( 



cos— - ^2(ga/9)cos6»J/(27ro 2 ), 



(44) 



<ltf(p,0,t)>=exp(- 



p 2 + (qa) 2 + 2pqa cos(<9 + 2irt/T c ) - 



( 



2irt A . N . \ . . o\ 
sin— + -^(qap)sm0)/(27ra ). 



(45) 



One can see from these expressions that at t = the x-component of spin density as a function of an angular variable 
has a maximum at 0q = tt (Fig. 4a). The maximum and minimum values of the function < T, y (p, 9, 0) > are in the 
vicinity of this point and differ in sign (see Fig. 5a). At t^ = Tr/A, as was shown above, we have two "classical" wave 
packets, Eq.(38), Fig. 2c. Correspondingly, 



1 



< X Xt y(p,0,t2) > = «(< E£ y (p,0,t2) > + < K$yMt2 +T cl /2) >), 



(46) 



that is illustrated by Fig. 4b and Fig. 5b. Note that such simple summation in Eq.(46) is only valid if the width of 
each of the wave packets is less than the cyclotron orbit length. 
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a 2 -£ y (/^t=0) 

0.00001 ^ ^ /77TTY 



(a) 




a 2 -L y (Aftt=T R /2) 
0.0006 



(c) 

FIG. 5: (Color online) The distribution of spin density < T, y (p,9,t) > for the wave packet, Eq.(24), with the parameters 
A/a = 0.1, qa = 5 at times: (a) t = 0, (b) t » T H /4, (c) £ » T^/2, (d) t « T R . 



IV. SPACE-TIME DYNAMICS OF THE MESOSCOPIC WAVE PACKETS FORMED BY THE 
POSITIVE AND NEGATIVE ENERGY EIGENSTATES 



In this Section we investigate the relativistic dynamics of the wave packet containing both the positive and negative 
energy states with \ k = 1. Correspondingly we restrict ourselves to the wave states with two nonzero components 
{^i, ^4}. As was noted in Section II the temporal dynamics of such states is governed by the part of Dirac Hamiltonian 
which is identical to a Jaynes-Cummings interaction. Besides of the general peculiarities of quantum dynamics of the 
system with the nonlinear energy spectrum (collapse and revival phenomenon), the existence of two energy bands 
leads to another effect - Zitterbewegung (ZB). This phenomenon, generating a highly oscillatory motion, is caused 
by the interference between positive and negative energy components of the wave packet. We will consider here the 
effect of ZB on the evolution of the electron velocity as well as on the spin polarization. 

Let at the time t — the initial wave packet is given by 



^(r,0) = ^ c (r) 



It> 



where 



Mr) 



1 / x 2 + y 2 . ixy\ 



is the wave function of coherent state, which can be written in the form 

^c(r) = / dp<p p (x)g(p) ^2c n+1 (j) n (y - y c ). 



(47) 



(48) 



(49) 



n=0 



The coefficients g(p) and c n in this expression are determined by Eqs.(22) and (23). Using Eqs.(20), (47) and (49), 
one can obtain 



^(r,r) = / dpip p (x)g(p)> 



£ 



71=1 



, dn<l> n - x {y - y c )\ t> \ exD /_ iw t) , b (K<t>n-i{y - Vc)\ t> I ,. T] 
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(50) 



"MnG/ " I/ c )| |> ; ^ V **"" ^ ^ V <W»(t/ " l/c)| l> 

Performing integration over p (see Appendix B, Eq.(B3)), we get finally 

^(r,r) = *i(r,r) (' ^) +* 2 (r,r) (j ^) , (51) 

where 

* l( ,,r) = ^(r)exp ( - f (^ - y + fa) ) £ ^g^f^ (cos^r) - 1*^I>), (52) 



n=l 



m 



,(r, r) = -^ c (r) exp ( - f (^ - „ + «)) £ ^(n -T i aV ^ " » + ^ (53) 



n=l 

As above, first of all we calculate the average velocity of the packet center 



vl(t) = c ^cFi^dr. (54) 



Substituting the expressions for the components of the wave function Eqs.(52), and (53) into Eq.(54) after the 
integration over the x, y coordinates we have 



v x (r) = c-exp( - (ga) 2 /2) V — — — ( cos(((/? n+ i - <p n )r) - cos((<p n+ i + ip n )r)\ 

a V / ^ 2 n i (n- l)!^ n y? n+ i \ / 



(55a) 



%M = c- exp ^ - (ga) 2 /2J ^ 2 n - 1 fn-l)! — ( sin (^ n + 1 ~ ^ n ) r ) ~ sin (On+i + ^n)r)j . (556) 

n=l ^ ' ^ n 

Note that these equations depend on not only on the energy difference <£ n +i — y? n as in Eq.(25), but on the sum 
(fn+i + ^ n also. The last dependence leads to the high frequency of the ZB oscillations. The dependence of the average 
velocity v^(r) on time for the packets parameters A/a = 0.5, qa = 10 (that corresponds to the dominant value of 
n = no — 50) is shown in Fig. 6. Like the case of the wave packet discussed in Section II, the behavior of this function 
is accompanied by the collapse and revival phenomenon. For the parameters of our wave packet the corresponding 
periods of cyclotron motion and revivals are T c \ w 126A/c and Tr « 26660A/c. In addition, the insert in Fig. 6 clearly 
demonstrates the ZB oscillations with peculiar frequency ujzb = 2(/? no c/A = 10.2c/A ~ 7.85 • 10 21 s -1 . Besides, these 
oscillations occur in short intervals near time moments tj~ = kT c i/2 (k = 1,2,.. .), that is connected with the feature 
of the space-time evolution of the wave packet probability density. At t > the initial wave packet (Fig. 7a) splits 
into two parts (Fig. 7b) with the amplitudes which are slightly different. These sub-packets containing the states of 
one sign of the energy rotate with cyclotron frequency in the opposite directions and meet each other every one-half 
of the cyclotron period (see Fig. 7c, 7d). The duration of ZB-oscillations is determined by the relation of the packets 
size to their relative velocity. Note that in the moment of time t = T c i/A the spin parts of the sub- wave packets with 
E n > and E n < differ slightly. Thus, the spatial part of the full wave function (50) is a linear superposition of the 
mesoscopic states (see Appendix C). At t ~ t^ the parts with positive and negative energy have significant overlap in 
the coordinate space. This is a necessary condition for the existence of Zitterbewegung (see insert in Fig. 6). 

At times t ~ T c \ the effect of the quadratic term in the Taylor expansion of the energy, Eq.(A2), is insignificant. 
For times much greater than the cyclotron period the dephasing appears in individual terms in Eq.(50) due to the 
term (n — no) 2 , that leads to the collapse of the sub- wave packets. At intermediate times T c \ <C t n < Tr, where 
t n w mTpi/n (m/n is irreducible fraction), the fractional revival of each sub-packets occurs. As a result, each sub- 
packet is decomposed into N = n(3 — (— l n ))/4 packets- fractions. Particularly, at t = Tr each of the two sub-packets 
(with positive or negative energy) restores at various points on the cyclotron orbit, that makes it impossible reshaping 
of the initial wave packet entirely. At time Tr/4 we should see four packets- fractions. However, this result is not 
consistent with the distribution of the probability density shown in Fig. 8a due to the significant influence of the term 
~ (n — no) 3 in the series (A2) for times t < Tr. This statement is illustrated in Fig. 8b, which shows at time t = Tr/4 
the probability density, taking into account the first three terms in the expansion (A2) only. 
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0.15 



■0.15 




FIG. 6: The dependence of the average velocity v x (t) versus time for the initial wave packet, Eq.(47), with the parameters 
A/a = 0.5, qa = 10. The oscillations of the average velocity v^(t) at times t > T c i are shown in the insert. 



Consider now the average value of the spin operator. Since the initial state of the wave packet, Eq.(47), belongs 



to one of the invariant subspaces (with A& = 1), the average spin projections S x 

S z = (h/tyfrz and Eqs.(51), (52), and (53) we obtain 



0. Using the definition 



3fr) = ? 



Em 



1 + 2n(\/a) 2 cos(2ip n T) 



<fl 



(56) 



where the coefficients c n are determined by Eq.(23). Corresponding dependence is shown in Fig. 9. It is clearly seen 
that the average value of the z-component of spin oscillates with the ZB-frequency ujzb ~ 2(/? no c/A and is accompanied 
by the phenomenon of collapse and revival. In this case unlike the previous definition, Eq.(30), the corresponding 
revival time is determined by the cyclotron period T rev = T c i/2&H In intervals between the ZB-oscillations the 
average value of S z (t) is not equal to zero. As mentioned above, the space-time evolution of the original wave packet 
(47) is described by the part of the Dirac Hamiltonian, Eq.(l) that is identical to the Jaynes-Cummings model in 
quantum optics. 7 So one can conclude that Eq.(56) is similar to the well-known expression for the population of the 
ground state of two- level atom exposed to the action of the classical electromagnetic field 1 . 



SUMMARY 



In this work we have analyzed the cyclotron dynamics of relativistic wave packets, moving in external constant 
magnetic field over cyclotron orbit. The components of electron wave function were calculated (see Eqs. (31), (52), 
and (53)), as well as, average velocities of the packet v^(t). It was shown that for the case when the wave packet 
consists of positive energy states only the effects of collapse and revival of electron and spin densities can be observed. 
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(a) 





(c) 

FIG. 7: (Color online) The distribution of the electron probability density for initial wave packet, Eq.(47), with the parameters 
A/a = 0.5, qa = 10: (a) t = 0, (b) t w T c ;/4, (c) t « T cZ /2, (d) £ w T cZ . 



|¥<;x.y.t=T R /4)l J 




|¥<x.y.t=T R /4) 




(a) 

FIG. 8: (Color online) Electron probability density for initial wave packet, Eq.(47), at qa = 10, A/a = 0.5, t ~ Tr/4: (a) for 
real energy spectrum, (b) for the spectrum obtained by taking into account the first three terms in the expansion (A2) only. 



On short-time t < T c \ (where T c \ is the cyclotron period) the coherent wave packet rotates with the classical 
frequency over the cyclotron orbit. We analyzed the complex evolution of spin densities and average spin components. 
The precession of spin densities one can see in Fig. 4, Fig. 5. We visualized N packets-fractions which form at times of 
fractional revivals t = mTji/n (N = n(3 — (— l) n )/4, as was shown by Averbukh and Perel'man 1 ). In our opinion, the 
dynamics of such structures can be studied experimentally by investigating the nature of electromagnetic multipole 
radiation and absorption. It was found that the effects of fractional revivals can not be observed via time-dependence 
of average velocity. 

We compared the dynamics of wave packet described above with the behavior of the wave packet containing the 
states with both the positive and negative energies (Sec. IV). In the last case the splitting of the initial wave packet 
takes place at t < T c \. These two sub-packets rotate with the cyclotron frequency in opposite directions and meet 
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A 9 - 



(t)/(f>/2) 



0.5 



-0.5 




FIG. 9: (Color online) The time dependence of average spin S z (in units of h/2) for initial wave packet, Eq.(47), for qa ■ 
A/a = 0.5. The oscillations of the average spin at short times t <C T c i are shown in the insert. 



10, 



each other every one-half of the cyclotron period that results in the high-frequency ZB oscillations. It was shown (see 
Appendix C) that of time t = T rev /2 = T c i/4 the electron spin and orbital degrees of freedom become disentangled. 
In this case spin states play the role of the internal states of the atom that has or has not radioactively decayed and 
the orbital states correspond to the cat in the original Schrodinger model 9 . At t — Tr each of the two sub-packets 
restores at different points of the cyclotron orbit. So that the full reshaping of initial wave packet is impossible unlike 
the case of wave packet including the positive energy states only. 
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Appendix A 



In this Appendix we find the approximate expression for the average electron velocity for the packet parameters 
qa ^> 1 and A/a <C 1 (we consider qa c± 5, A/a ~ 0.1). Corresponding wave function (20) represents the quantum 
superposition of the states sharply peaked around large enough central value n = uq. In fact, for (ga) 2 /2 ^> 1 the 
Poisson distribution \c n \ 2 (see Eq.(23)) can be replaced by Gaussian one 



<qa) 2 /2 



m n 



V27rn 



exp 



(n -n ) s 
2n 



(Al) 
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where no — (qa) 2 /2 c± 12 and dispersion An ~ v^q. So, we expand the dimensionless energy spectrum cp n as a power 
series in (n — no), restricting ourselves to the quadratic terms^ 

¥n = (Pno + ^n ( n ~ n o) + -^(jl ~ n ) 2 . (A2) 

For the given values of the packet parameters (p n is close to 1 and the coefficients before the exponents in Eq.(25) can 
be estimated, in good approximation, as 



y 2(n + l)y? n+ i \/2 ' V ^n 
that allows us to carry out the sum over n in Eq.(25). The result is 

*Mt) = — exp I - lqasm—^—J Jx 



^ +1 -1 _,A/a MT + T^^ (A3) 



cos I ^ I cos ( <,r - ^ o ( -^ - 1 ) T + -^- sin^ r I , (A4) 



„ / (ggr ,\ , (gq) • » \ 

vZ{t) = — exp I - (qai 
m \ V 



, <n 2 



r,sm— ^-] \y. 

COS 



v~2~ J sm v^ n ° r " ^ no v^T~ " ) T + ~l~~ sm ^ n ° r J • ( ^ 

As clearly seen from these expressions the velocity oscillations reduce exponentially with time. The corresponding 
damping time is determined by Eq.(29). 

Appendix B 

i) Below we find more detailed expression for wave function (24). In order to do that we perform the integration 
over p. The corresponding integral has the form 

+oo 

Qk= <p P (x)g(p)(l>k(y-yc(p))dp, (Bi) 

— oo 

where (f p (x) = (l/y/27rh) ex.-p(ipx/h), where g(p) is determined by Eq.(22) and 4>k{y — Vc(p)) is the well-known linear 
harmonic oscillator function. Putting these expressions into Eq.(Bl) and using Gaussian transformation for Hermite 
polynomials Hk(y)r^ 

+oo 

^ j e-^-y^ H k {y) dy = (2x) k , (B2) 



we finally obtain 



v^ 



((y - Q a2 — ix)/a) k /2ix(y + qa 2 ) — x 2 — (y — qa 2 ) 2 ' 



^2 k + 1 k\ira 2 



/Zix(y + qa") - x" - \y-qary\ , . 

{ ^5 J ' ( 53 ) 



Substituting the above expression into Eq.(24) and using polar coordinates (32) one gets Eq.(31). 

ii) Let us now calculate the electron wave function ip(r, t) = ipcl( r > t) at the short times t <C Tjj when we can neglect 
the quadratic term in the expansion (A2). Then recollecting that for a given wave packet the coefficients d n and b n 
can be approximated by the expressions d n ~ 1 and b n ~ y/n/2A/a, we find after summation from Eq.(31) 



17 

where 

S(u) = exp ( — (cos u — i sin u) J (B5) 

and M{p,8) is determined by Eq.(33). The solution (B4) describes the coherent weakly-relativistic wave packet 
propagating along the cyclotron orbit without changing its shape. In particular, at a = /3 the corresponding probability 
density has the Gaussian form 

\* d {p,0,t)\* = l + ^ {p2 t {Qa)2) exp (- p2 + ^ + 2pqaC °< 9 + 2 ^ T -\ (B6) 

2ira z V 2 / 

It is clearly seen that the center of the wave packet corresponding to the angle 

rotates with cyclotron frequency cu c = 2ir/T c similar to the motion of a classical electron placed in a magnetic field. 

Appendix C 

Below we show that the wave packet discussed in Sec. IV (Eqs.(47), (48)) at to = T c i/4 is a coherent superposition 
of mesoscopic states&i^. To do it we consider the time evolution of the initial state 

V>+(r,0) = Vc(r)|+>, (CI) 

where the "orbital" wave function of the coherent state ip c ( r ) is given by Eq.(48) and the initial Dirac spinor |+ > 
has the form 

l+>=<('| ) > )+^o(|] ) > ) (C2) 

with the coefficients d no and b no are given by Eq.(14). Here as above we assume that the mean value of n in Eq.(49) 
is large enough: n = (qa) 2 /2 = no ^> 1. 
The state (CI) evolves as 



^+(r,r) = / dpcp p (x)g(p) ^ c+ I 

J Ln=l ^ 

where the coefficients cj are 



dn^n-liV ~ Vc)\ t> \ e -zcp n r 



Wn-l(|/-|/c)| t> ] e %ip n r 

x d n (f) n (y - y c )\ i> 

n=0 



b„<Pn{y-yc)\ l> 

(C3) 



c~ = c n d no b n + c n+ ib no b n . (C4) 

Using the weak dependence on n for the coefficients d n and 6 n and also using the expression for c n , Eq.(23), we obtain 
from (C4) at n ~ no ^> 1 

c+ = c n , c~ = 0. (C5) 

Then it is not difficult to show that 

iMr.r) = (V^o^„ ('],>)+&»„ ( | °>) )^c+(r,r), (C6) 



where 



^c+(r,r)= dp(f p (x)g(p)^2c n+1 (t) n (y - y c )e lifnT . (C7) 

^ n=0 
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In the same way one can compute the evolution of the initial state 

^_(r,0)=Vc(r)|->, (C8) 

\->= b no( l l > )-dno( l l > Y (C9) 



V_(r,r) = ( e^o^o ( ' ^ ) - d no ( , ^ ) ^-(r.r), (CIO) 



*>)-<* ( ° 

o J ™° ^ \i> 

where ip c -( r , t) is determined by Eq.(C7) in which exp(— itp n r) should be replaced by exp(iy> n -r). Note that the initial 
spin functions are orthogonal: < +|— >= 0. 

Thus, the initial wave packet, Eq.(47), can be represented as a superposition 

Mr) ( l l > )= d no ^+{v, 0) + &„ ^-(r, 0), (Cll) 



so that 

^(r,r) =^ no ^+(r,r)+6 no ^-(r,r) (C12) 

describing two sub-packets which at t < T c \ = (27r/cp f no )X/c rotate counterclockwise (^+( r ? r )) an d clockwise (^-(r, r)) 
in the xy— plane. At the moment of time to = T c i/A as follows from Eqs.(C6) and (CIO) 



< + -> 



t=t 



2n (A/a) 2 
+ 2n (A/a) 2 



0.98 (CIS) 



for no = 50 and A/a = 0.5, and we conclude that the spin parts of wave functions ^±( r ? t) differ slightly. That means 
that "orbital" part of full wave function (C12) at that time is a linear superposition of the mesoscopic states ^ c +(r, to) 
and i/j c -(r,to) (the so-called "Schrodinger cat"). The corresponding probability density |^(r,£o = T c i/4)\ 2 is shown 
in Fig. 7b. The generation of coherent superposition of mesoscopic states for the Dirac particle placed in the magnetic 
field was predicted in Ref. 6 . 
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